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The "^Cu Knight shifts and g factors for the normal state of YBa2Cu307 in tetragonal phase are theoretically 
studied in a uniform way from the high (fourth-) order perturbation formulas of these parameters for a 3d:' ion 
under tetragonally elongated octahedra. The calculations are quantitatively correlated with the local structure 
of the Cu^+(2) site in YBa2Cu307. The theoretical results show good agreement with the observed values, 
and the improvements are achieved by adopting fewer adjustable parameters as compared to the previous 
works. It is found that the significant anisotropy of the Knight shifts is mainly attributed to the anisotropy of the 
g factors due to the orbital interactions. 
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1. Introduction 

YBa2Cu307 has become an important topic in the microwave magnetism i2i] as well as 
nuclear quadrupole resonance (NQR) and nuclear magnetic resonance (NMR) 0, [^l researches. In 
particular, the Knight shifts as the results of NMR experiments can yield important information 
about the local environments and electronic distribution around the magnetic nuclei. For example, 
the Knight shifts and K±^ as well as the anisotropic g factors and g± were measured 0, 0] 
for the tetragonal Cu^^(2) site in YBa2Cu307 of normal state at room temperature. However, the 
theoretical studies [6j] on the Knight shifts seem unsatisfactory. Firstly, the previous calculations 
of the Knight shifts were generally based on the simple second-order perturbation formulas, while 
the contributions from the higher (third and fourth) order perturbation terms were not taken into 
account. Secondly, the analysis on the Knight shifts was not correlated to the local structure of 
the magnetic site, but was treated by introducing various adjustable parameters (e.g., the related 
energy separations). Finally, the g factors were not quantitatively treated in a uniform way. In 
order to study the Knight shifts and the g factors of YBa2Cu307 to a better extent, the high 
(fourth-) order perturbation formulas of the Knight shifts and the g factors for a ion located 
in tetragonally elongated octahedra are applied to the tetragonal Cu^^(2) site of YBa2Cu307 in 
this work, and the local structure of this site is quantitatively involved in the calculations. 

2. Calculation 

The anisotropic Knight shifts and g factors may be ascribed to the tetragonal Cu^^(2) site in 
YBa2Cu307, which is coordinated to five oxygen ions forming a tetragonally elongated octahedron 
(i.e., one of the apical ligands in the oxygen octahedron is moved to infinity) [8j] As regards a 
Cu'^~^{3dP) ion in tetragonally elongated octahedra, its original cubic '^Eg ground state would 
be separated into two orbital singlets ^Big (or jx^ — y^)) and ^Aig (or \z^)), and the former is 
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the lowest level. Meanwhile, the original cubic '^T2g exited state may spilt into an orbital singlet 
^B2g (or \xy)) and a doublet '^Eg (or \yz) and \xz)) [9]. In order to make uniform investigations 
of the Knight shifts and the g factors as well as to overcome the shortcomings of the previous 
study ^ based on the simple second-order perturbation formulas, the high order perturbation g 
formulas [lOj including the third- and fourth- order perturbation terms for a tetragonally elongated 
octahedral 3(P cluster are adopted here. Applying the Macfarlane's perturbation loop method [ll|, 
the perturbation Hamiltonian for an orthorhombic 3(P cluster under external magnetic field can 
be written as follows: 



H' 



(1) 



where Hgo, Hzc and iJhf are, respectively, the spin-orbit coupling interactions, the Zeeman term 
and the hyperfine interactions, with the corresponding spin-orbit coupling coefficient ^, the orbital 
reduction factor k and the dipolar hyperfine structure parameter P for a ion in crystals. Uti- 
lizing the perturbation method the perturbation formulas for an orthorhombically elongated 
3(P cluster can be derived as follows lldl: 
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here gs{~ 2.0023) is the pure spin value, k is the orbital reduction factor. C, is the spin-orbit coupling 
coefficient for the Sd^ ion in crystals, which can be expressed in terms of the corresponding free- 
ion value Co as « fcCo- It is noted that when the third- and fourth -order perturbation terms 
are neglected, i.e., only the first and second terms in the right sides of equation ([T]) are reserved, 
the above formulas may be reduced to those of the previous work Ei and E2 are the energy 
separations between the excited ^i?2g and ^i?g and the ground ^-Big states [12l |: 
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Here Z?q is the cubic field parameter, and and Dt are the tetragonal field parameters. The 
local structure of the five- fold coordinated Cu^~''(2) site in YBa2Cu307 can be described as the 

2.426 A) along c axis and the four perpendicular distances 
i|. Thus, the tetragonal field parameters can be determined 



parallel Cu — O distance R\\ (« 
R±i~ 1.939 A) along a and b axes 
from the superposition model (l2| : 
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Here ^2 ~ 3 and t4 « 5 are the power-law exponents [12|. A2{Ro) and A4{Ro) are the intrinsic 
parameters, with the reference bond length Rq taken as the average Cu^"*" — O^^ distance, i.e., 
Rq = R = 4- 4i?j_)/5. For octahedral Sd" clusters, the relationships ^4(i?o) ~ (3/4)-Dq and 
A2{Ro) ~ 10.8A4(i?o) are proved valid for many systems and reasonably applied here. 

According to the optical spectra for Cu^"'' in oxides the cubic field parameter Dq « 
1260 cm~^ and the orbital reduction factor k « 0.76 can be obtained for the studied system by 
fitting the observed d-d transitions. The spin-orbit coupling coefficient is determined using the free- 
ion value Co(~ 829 cm"-'^ [14]) for Cu^^. Substituting these values into equation ([T]), the theoretical 
g factors are calculated and shown in table 1. 
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Now we turn to the investigations of the Knight shifts. From the relationships between the 
g factors and the Knight shifts and the perturbation procedure similar to that in Ref. the 
high order perturbation formulas (i.e., those containing the third- and fourth- order perturbation 
contributions) of the Knight shifts for the tetragonal Cu^~''(2) site in YBa2Cu307 can be expressed 
as: 
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In the above expressions, x — ■^)3d-^AMB- Here {r^^)3di~ 4.68 a.u. [7]) is the expectation 
value of inverse cube of the 3d radial wave function of Cu^^ in the system studied. A'a is the 
Avogadro's number, /ie is the Bohr magneton. Similarly, since the third- and fourth-order pertur- 
bation terms are neglected, i.e., only the first terms in the right sides of equation ([SJ are reserved, 
the above formulas are reduced to those of the previous work Substituting these values into 
equation (O, the Knight shifts are calculated and collected in table 1. For comparisons, the the- 
oretical results of the previous work Q based on the simple second-order g factors and various 
adjustable energy separations are listed in table 1. 



Table 1. The Knight shifts and the g factors for the Cu'^+(2) site in YBa2Cu307 of normal state. 
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a Calculations based on the simple second-order g factors and various adjustable energy separations 
in the previous work Q. 

b Calculations based on the uniform high order perturbation formulas and the local structure of the 
tetragonal Cu'^'''(2) site in YBa2Cu307 in this work. 



3. Discussion 

From table 1, one can find that the theoretical Knight shifts and g factors of this work show 
reasonable agreement with the experimental data, and these parameters are also suitably explained 
in a uniform way. 

In the present calculations, the relationships between the Knight shifts (and the g factors) and 
the local structure of the tetragonal Cu^''" site in YBa2Cu307 are quantitatively established from 
the superposition model, and the shortcoming of the previous work [6| based on various adjustable 
energy separations is therefore overcome. Meanwhile, the high order perturbation formulas con- 
taining the third- and fourth-order perturbation contributions in this work are also superior to the 
simple second-order ones in the previous study. It is noted that a crude estimation of the contribu- 
tions of the still higher (fifth) order perturbation terms yields ~ C^/(i?^i?2), in the order of 10^^ 
and being safely negligible. In addition, the previous results (Cal.") comparable with the observed 
values in view of the experimental uncertainties may be illustrated by the fact that the calculation 
errors happened to be canceled by the various adjustable energy separations 
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The positive anisotropy — K± is consistent with the positive anisotropy g\\ — g±. This 
can be ascribed to the approximate Hnear relationships between the g factors and the Knight 
shifts. Further, the above anisotropics may be attributable to the local tetragonal elongation of the 
Cu^^(2) site, i.e., the relatively longer parallel Cu-0 bond as compared with the four perpendicular 
ones. Thus, the ground ^-Big state and hence the anisotropic orbital angular momentum and spin 
interaction between electron and nuclear can be understood. The calculated anisotropy — 
K±{Ki 1.081 in Cal.'') of this work is somewhat larger than that (« 0.996 in Cal.°) of the previous 
work |6|. This can be illustrated by the fact that the present anisotropy g\\ — g±{^ 0.251 in Cal.'') 
is smaller than the previous result (« 0.27 in Cal.'^) [6] and hence slightly worse as compared with 
the observed value (w 0.30 (20) 0,01). The above discrepancies are ascribed to the approximations 
of the theoretical model and formulas as well as to the calculation errors in this work. In view of 
the experimental uncertainties of the observed A'n — K±^ and — g±, the present calculations can 
still be regarded as suitable. 
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TeopeTMHHe BUBHeHHJi scysiB HawTa ^^Cu 
B ocHOBHOMy CTahi YBa2Cu307 

M.K. KyaHrS, C.-M. B^, S.-V. yKan^, B.T. CoHrS 

0aKy/ibTeT npwKJiaflHoT cfjisnKH, VHiBepci/iTeT e/ieKxpoHiKi/i i jexHo^iorlCi Ki/iraio, 610054 HeHTfly, Ki/iraii 
Mi>KHapoflHMi^ L4eHTp 4>i3WKM MaTepiajliB, Ki/iTaWcbKa aKa^eMl;; HayK, 110016 LUeHbaHr, KwraM 



3cyBM Hai^Ta ''^Cu i g 4>aKTopw b ocHOBHOMy cxaHi YBa2Cu307 b xexparoHajibHii^ dpa3\ BUBHaKiTbCJi xeope- 
TWHHO B eflWHOMy niflxofli 3 4>opMyji xeopiT sBypeHb Bi/icoKoro (neTBepToro) nopaflKy liwx napaMexpiB p,r\9f Scfi 
ioHa B TerparoHa/ibHO BUflOBMeHowy OKxaeflpi. PospaxyHKi/i KiiibKicHO KopejiKDKDTb 3 /lOKa/ibHOKD crpyKTypoKD 
Cu^'^{2) By3Jia b YBa2Cu307. TeopernHHi pe3y/ibTaTM fleMOHcrpyHDTb flo5pe y3rofl>KeHH5i 3i cnocTepe>Ky- 

BaHMMH 3HaHeHH!RMW, i HOKpa LU,eHH51 flocjiraeTbcs LUJiaXOM BI/lKOpi/lCTaHHa MeHUJoT Ki/lbKOCri niflrOHOHHI/IX 

napaMexpiB Hi>K y nonepeflHix po6oTax. SHaCifleHO, L140 SHaHHa aHisorponin 3cyBiB Hawra e b ocHOBHOMy 
aTpw6yTOM aHiBorponif g cJjaKTopiB 3aBA^Ki/i B3aeMOfliT. 

K/iiOHOBi c/ioea: e/ieKTpoHHuH napaMarneTi/iHHi/i^ pe30H3HC, affepni/i^ MarniTHi^H pesoHaHC, 3cyB HaMTa, 
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